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FORMULATION OF IMAGE QUALITY PREDICTION CRITERIA 
FOR THE VIKING LANDER CAMERA 

By Friedrich O. Huck, Daniel J. Jobson, Edward J. Taylor, 
and Stephen D. Wall 
Langley Research Center 

SUMMARY 

Image quality criteria are defined and mathematically formulated for the prediction 
computer program which is to be developed for the Viking lander imaging experiment. 
The general objective of broad-band (black and white) imagery to resolve small spatial 
details and slopes is formulated as the detectability of a right-circular cone with surface 
properties of the surrounding terrain. The general objective of narrow-band (color and 
near -infrared) imagery to observe spectral characteristics is formulated as the mini- 
mum detectable albedo variation. The general goal to encompass, but not exceed, the 
range of the scene radiance distribution within a single, commandable, camera dynamic 
range setting is also considered. 


INTRODUCTION 

One of the primary experiments of the two Viking lander missions to Mars in 1975 
is imagery, two facsimile cameras being used on each lander. The primary purpose of 
the imaging experiment is to spatially and (to a lesser extent) spectrally characterize 
the terrain which surrounds the landers . Optimization of this experiment depends not 
only on the camera design but also on the imaging strategy. It is, for example, generally 
known from the Lunar Orbiter and Surveyor missions to the Moon and the Mariner mis- 
sions to Mars that image contrast of surface detail is critically dependent on illumination 
and viewing geometry. 

Preflight predictions of image quality have been made by Rindfleisch and Willingham 
(refs. 1 and 2) for the Ranger missions, and by Keene (ref. 3) and Huck (ref. 4) for the 
Lunar Orbiter missions. These investigations have shown that a careful formulation and 
evaluation of the interaction between camera characteristics, mission constraints, and 
known or estimated surface properties can contribute significantly to a successful 
imaging strategy. Although these investigations and the experiences gained from the 
Surveyor imaging experiments are helpful in establishing an imaging strategy for the 



Viking lander experiment, an image prediction analysis is nevertheless needed. This 
need arises because the performance characteristics of the facsimile cameras differ 
from previously used television and film cameras, and because mission constraints differ 
from previous space imaging experiments. 

Of primary concern to both the Viking lander camera design and imaging strategy - 
as should prove typical for exploratory planetary imaging experiments in general - is 
the unknown nature of surface properties at the scale of observation and the limited data 
transmission capacity due to the long distances and spacecraft weight and power con- 
straints. The camera must, therefore, have a good sensitivity over a wide dynamic 
range whereas the number of encoding levels of surface radiance must be kept low. To 
accommodate these two conflicting requirements, the Viking facsimile camera uses 6-bit 
encoding with 192 partial dynamic ranges which can be selected by 6 gains and 32 offsets, 
These partial dynamic ranges have a linear transfer function to provide images with good 
radiometric accuracy, rather than a logarithmic transfer function which would insure 
superficially pleasing image quality. 

In addition, the Viking lander cameras feature four imaging modes (survey, high 
resolution, color, and near infrared) and two scan rates (matched to the data transmis- 
sion rates of 16 000 bits per second to the Orbiter and 250 bits per second directly to 
Earth); these features result in various performance characteristics and require differ- 
ent viewing geometries. The cameras also have a single-line scan mode which permits 
an increase in radiometric sensitivity by repeatedly scanning a single line in the scene. 

Preflight predictions of image quality should be conducted both experimentally and 
analytically. The primary advantage of experimental predictions is that results that 
appear in the form of images will principally be used for interpreting the results of the 
actual investigation, whereas the primary advantage of analytical predictions is that a 
wide range of conditions can readily be simulated on the computer. Thus, a handbook of 
images obtained under various imaging conditions and with various camera settings can 
probably serve the experimenter best if anticipated situations (for example, of Mars 
optical surface properties and camera performance characteristics) are encountered, 
and a flexible software program can probably serve him best if unanticipated situations 
are encountered that require a rapid response. The pursuit of both approaches not only 
combines these advantages but also provides a means for checking the results of one 
approach with the other. 

The purpose of this paper is to provide the necessary analytical basis for devel- 
oping the image quality prediction computer program for the Viking lander imaging 
experiment. Presented are definitions and mathematical formulations of pertinent image 
quality criteria as functions of optical surface properties, illumination and viewing geom- 
etry, and camera characteristics. It is important to recognize that an image quality 
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prediction computer program based on this analysis must be carefully compared with 
experimental preflight investigations before becoming operational. 

SYMBOLS 

o 

A area of a cone section, meters 

b normalized diameter of approximated circular area of a cone section 

(see eqs. (16)) 

c normalized diameter of cone base, 2R C/ //3L 

D diameter, meters 

f lens focal length, meters 

g phase angle (see fig. 2), degrees 

H height, meters 

I current, amperes 

i integer 

J n ( ) nth order Bessel function where n = 0, 1, 2, . . . 

k spatial frequency, line pairs per meter 

L distance from camera lens in object space (see fig. 5(b)), meters 

L(u,s) lens spatial frequency response 

t distance from camera lens in image space (see fig. 5(b)), meters 

m number of binary digits used for encoding shades of gray, bits 

spectral radiance, watts/meter^-micrometer-steradian 
N unit vector normal to surface 
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n 

p ( v a s ) 

p x 

p x( bv a s ) 

R 

% 

Sx 

S/N 

s 

u 

V 

a 

P 


number of overlapping line scans 

spatial frequency response of circular photosensor aperture 
spectral radiant power, watts /micrometer 
spatial frequency distribution of circular area of a cone section 
radius, meters 

photosensor responsivity, amperes/watt 

solar irradiance above Martian atmosphere, watts/meter^ -micrometer 
ratio of average signal to root-mean-square noise 
dimensionless variable for spatial frequency k (see eq. (15)) 
dimensionless variable for defocus (see eq. (20)) 
dimensionless variable for radius (see eq. (15)) 
target slope, degrees or radians 

instantaneous field of view or angular resolution (see fig. 5(b)), degrees or 
radians 



6 azimuth cone angle increment (see fig. 3), degrees or radians 

e angle between emitted radiation and normal to surface (see fig. 2), degrees 

or radians 

£ azimuth angle between object slope and incident radiation (see fig. 2), 

degrees or radians 


9 azimuth angle between incident and emitted radiation (see fig. 2), degrees or 

radians 
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l angle between incident radiation and normal to surface (see fig. 2), degrees 

or radians 

X wavelength, meters 


* ■ |(?sh + f-sh - s ) 


n(x) top hat function ^=1 if |x| s I; =0 elsewhere^ 

spectral reflectivity of surface (normal albedo) 
cr population deviation 

T x( t o) spectral transmissivity of atmosphere 

g spectral transmissivity of optics 

4> illumination scattering function 

4/ azimuth angle between object slope and emitted radiation (see fig. 2), degrees 

or radians 

angle defining numerical aperture (see fig. 5(b)), degrees 
a; angle, radians 

Subscripts: 

a photosensor aperture 

c cone 

e electronic 

i integer 

l lens 
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n noise 

o flat surface 

q quantization 

s signal 

sh shadow 

v visible 

+ brighter than flat surface 

darker than flat surface 

The symbol A in front of a parameter indicates a differential of that parameter. 
The bracket ( ) around a parameter indicates a criterion which is to be estimated. 
Bar over a symbol denotes an average value. A circumflex (~) over a symbol denotes 
a unit vector. 


DEFINITION OF IMAGE QUALITY CRITERIA 

The Viking lander camera (see ref. 5 and the appendix for details) has four imaging 
modes: survey, high resolution, color, and near infrared (IR). In the broad-band survey 
and high -resolution imaging modes, it is generally desirable to record small spatial 
details and slope variations. In the narrow -band color and IR imaging modes, it is gen- 
erally desirable to record spectral variations. In addition, it is generally also desirable 
in all four imaging modes to encompass, but not exceed, the complete range of radiance 
variations in the scene with a single dynamic gain setting. The purpose of this section is 
to define image quality criteria which will permit the mathematical formulation of these 
three objectives: 

(1) Knowledge of the statistical distribution (histogram) of surface radiance would 
allow the selection of an optimum camera dynamic range setting (that is, gain and offset 
as illustrated in the appendix) for the complete or a particular part of the radiance range. 
But the actual radiance distribution will, of course, not be known until after video data 
have been received from the lander. Until such data have been received, it is proposed 
that a minimum, mean, and maximum value of the radiance distribution be estimated, the 
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albedo and illumination scattering function of the scene being assumed to be uniform over 
the landing site. Clearly, the minimum surface radiance occurs in shadows. Since 
atmospheric scattering on Mars is small within the spectral range of the camera silicon 
photosensors (excepting, of course, a dust storm), shadows will exhibit very little 
radiance - certainly much less than the smallest commandable camera offset. (See 
appendix.) Hence, for the purpose of selecting a camera offset, minimum radiance may 
be defined as zero. The mean radiance is defined as the radiance of a level surface for 
a given viewing geometry. The maximum radiance is defined as the radiance of' a surface 
area with a slope and slope orientation which produces highest surface reflectance, again 
for a given viewing geometry. Hence, mean and maximum radiance estimates are func- 
tions of surface albedo and illumination scattering function. 

(2) The capability to resolve spectral variations is defined as the minimum detect- 
able albedo variation. This minimum detectable variation is taken to be that difference 
in albedo which results in a ratio of signal to root-mean-square noise of 3, a level sur- 
face being assumed. 

(3) The capability to resolve small spatial details and slopes is defined here as the 
minimum detectable cone diameter and cone slope with respect to a level surface. This 
criterion is based in part on a figure of merit for image quality which was first proposed 
in reference 1. 

A right-circular cone with surface properties of the surrounding terrain seems 
intuitively representative of many features and has no preferred surface orientation azi- 
muthally about its axis. If the cone angle is chosen to be steep, a condition yielding high 
surface contrast, then the detectability of this target becomes primarily a measure of 
the camera capability to resolve small detail. If the cone angle is chosen to be shallow, 
a condition yielding low surface contrast, then the detectability of the same target 
becomes primarily a measure of the camera capability to resolve small slopes. 

The cone is probably also the simplest shape for this application. But it would, 
nevertheless, be unnecessarily complex to translate rigorously cone-surface-radiance 
variations into image grey-scale variations in order to estimate the detectability of a 
cone. Instead, the average radiance is calculated separately over those regions of the 
cone which are brighter than a level surface background, and those regions, including a 
shadow (if present), which are darker than the background. The two average reflectance 
values and their corresponding areas are then used as approximate target characteristics 
in order to formulate an expression for the ratio of signal to noise for the image of a 
cone. 
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MATHEMATICAL FORMULATION OF IMAGE QUALITY CRITERIA 


Figure 1 outlines the analytical approach taken to formulate the image quality cri- 
teria which have just been discussed. Each block of this flow diagram represents a sec- 
tion of the analysis. 

The first section specifies the angular relationships necessary to account for the 
dependence of surface reflectance on surface slopes and illumination and viewing geom- 
etry. The next section presents the geometry of an upright cone, which was selected to 
represent small spatial detail and slope variations. Reflectance variations of the cone 
are approximated in the third section by an average reflectance value for those sections 
of the cone which have a higher reflectance than a level background surface and another 
average reflectance value for those sections which have a lower reflectance (including 
shadow, if present). The fourth section formulates the resultant spectral surface radi- 
ance and the radiant power incident on the camera, and the fifth section formulates the 
camera response to this input signal. The sixth section describes the noise generated 
by the photosensor and preamplifier, and the quantization error caused by encoding the 
electrical signal for digital transmission. The last section formulates the desired image 
quality parameters. 


Viewing Geometry 

Figure 2 defines the angular relationships between surface slopes and illumination 
and viewing geometry. The illumination scattering function <^(e,i,g), which accounts 
for the dependence of surface reflectance on this geometry, is a function of the angle i 
between incident radiation and surface normal, the angle e between emitted radiation 
and surface normal, and the phase angle g between incident and emitted radiation. 
These three angles can be determined from the following relationships: 

cos g = cos e 0 cos i 0 + sin e 0 sin l q cos 9 (1) 

cos i = cos a cos i 0 - sin a sin i 0 cos £ (2) 

cos e = cos a cos e 0 - sin a sin e 0 cos \(/ (3) 

where e Q and i Q are the emission and incident angles, respectively, for a level sur- 
face (with normal N 0 ), a is the angle of a surface element with respect to a level 
surface (and with normal n), 6 is the azimuth angle between incident radiation and 

camera optical axis, and £ and 4/ are azimuth angles of the plane formed by the 
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normals N and N 0 with respect to incident radiation and camera, respectively. The 
three azimuth angles are related by 

6 + 5 + ify = 2tt (4) 


Target Geometry 

Geometric relationships required to describe the reflectance variations of a right 
circular cone have been derived in reference 1 and are generalized here to include 
viewing geometries for that part of the cone, including its shadow (if present), which is 
obscured from the camera view by the cone itself. Only an upright cone is considered 
here, whereas reference 1 also considers an inverted cone. 

Figure 3 illustrates the geometry of an upright cone with a base of radius R c , a 
height H, and a slope a so that 


= tan_1 (° < “ < f) ( 5 ) 

An area element about a generator of the cone and with normal N is given in terms of 
the azimuth angle increment 5 as 


AA = 


6R C 2 

2 cos a 


( 6 ) 


If t 0 = 90° - a, a shadow is cast by the cone on a level surface (as shown shaded 
in fig. 3(a)) for -£ s h < ? < where Cg^ = cos"^(l/tan a tan i 0 ); this shadow has an 
area 

•^sh = Rc (t an ?sh " ^sh) 

Similarly, if e Q g 90° - a, a "visual" shadow is formed (fig. 3(b)) which obscures 
part of the cone and cone shadow (if there is one) from the camera. It is assumed that 
the rays of light intercepted by the camera are parallel instead of convergent since only 
small targets are of interest and the convergence angle is therefore small. Under this 
condition, area elements AA of the cone not visible to the camera are defined by the 
range where = cos'^(l/tan a tan e 0 ). The shadow area still 

visible to the camera becomes (see fig. 4) 
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(7b) 


Ay,sh _ Rc^(t an ?sh " ?sh “ t- 311 £ + £) 

where * = |-p sh + «^ sh - 0). 

Several geometries must be accounted for when equation (7b) is used. The azimuth 
angle 8 between incident radiation and camera optical axis must always be chosen so 
that 6 < 77. When 6 ^ the "visual" shadow does not overlap any part of the 

illumination shadow, and equation (7b) reduces to equation (7a). When 0 + £ sh g i// gh , the 
illumination shadow is completely contained within the visual shadow, and the angle £ 
must be set equal to £ sh so that A sh = 0. When 6 + i^ stl S the visual shadow is 
completely contained within the illumination shadow, and the angle £ must be set equal 
to rf/ sh . 


Approximate Target Characteristics 

Target and background are considered to be of the same material so that the varia- 
tion in cone reflectance with azimuth angle is dependent only on the illumination scatter- 
ing function. But it would still be exceedingly difficult to translate this reflectance vari- 
ation exactly into an image signal. The reflectance variation is therefore approximated 
by an average reflectance of that part of the cone which has a higher than background 
reflectance and another part which has a lower than background (including shadow) 
reflectance. 

The total projected area of all cone sections which have a higher than background 
reflectance and are visible to the camera may be expressed as 

A + = AA £ cos Cj = £ cos €j (*!<€, i,g) > *o) (8) 

i i 

where (p Q = <t>(^ 0 , L Q ,g )• The average value of the illumination scattering function of this 
area becomes 

^ 0 i ( e , t ,g)cos e i 

0 + = ^ ^ 0i(e,i,g)cos e i = -i — — (<p^(e,L,g) > <f> 0 ) (9) 

+ i 2 C0Sei 

i 

Similarly, the total projected area of all cone sections, including the shadow cast 
by the cone on a level surface (if present), which have a lower than background reflec- 
tance and are visible to the camera may be expressed as 
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( 10 ) 


A_ = AA cos + A s ^ cos e Q 
i 

rars I cos £ i + c sh - ?sh - ““ « + «) cos e o 

i 

(0i(€,i,g) < 0 O ) 

The average value of the illumination scattering function of this area becomes 



4 >_ = ^ ^ 0i(e,<-,g)cos €. 
i 

6^ <^ i (e,(.,g)cos e i 


6 ^ cos €^2 cos a cos e Q (tan £ sh - ? gh - tan ^ 
i 

(^(e,t,g) < 4> 0 ) 

Reflectance from the shadow is assumed to be negligible. 


(ID 


Spectral Radiance and Radiant Power 
The spectral radiance of a surface is generally given by (ref. 6) 


N x = fVx( l o)<W £ -‘-f:> < 12 > 

where is the solar irradiance, T x{ L o) * s t ransm i ss i v ity of the atmosphere, and 
p^ is the reflectivity (normal albedo) of the surface, (if p^ is defined as bond albedo, 
then the factor - should be replaced by (See ref. 7.)) 

7T 07T / 

The spectral radiant power reaching the camera photosensor is (ref. 8) 


p x * (DVd/tx s Nx = j| ^V s x T x( l o) I’xXstx^-s') 


(13) 


where /3 is the camera instantaneous field of view or angular resolution; D^, the 
objective lens diameter; and g , the transmissivity of the optics and spectral filters. 

To obtain the spectral radiant power Px o °1 a level surface, the function 
<p^(e,L, g) in equation (1 3) may be replaced by (f> Q . To obtain the spectral -radiant-power 
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difference between the level surface and the brighter than background cone area AP^ , 
and the darker than background area AP A the function <^(e,L, g) may be replaced in 
equation (13) by A0 + = 4> + ~ 4> Q and A$_ = <p Q - 4>_, respectively. 

Imaging Process 

The next step is to formulate the method by which the facsimile camera translates 
the radiance variations formulated in the foregoing sections into an electrical signal. A 
basic configuration of the camera is shown in figure 5. Radiation from the scene is 
captured by the scanning mirror and objective lens and projected onto a plane which con- 
tains the photosensor aperture. The photosensor - in this case, a silicon photodiode - 
converts the radiation falling on the aperture into an electrical current. As the mirror 
rotates, the imaged scene moves past the aperture and permits the aperture to scan ver- 
tical strips. The camera rotates slowly in azimuth so that the entire scene of interest is 
scanned. 

The line-scan imaging process of the facsimile camera has been formulated in ref- 
erence 8, and all pertinent optical characteristics have been presented in reference 9. 
Results from these references are used here in the notation of the present report without 
detailed derivations, and effects of line-scan sampling, such as possible aliasing (ref. 8), 
are not considered. 

The spectral radiant power P x reaching the photosensor aperture from a sur- 
face of uniform radiance is converted by the photosensor into a current of magnitude 

’o = r p x, 0 R x ■» < h > 

where R^ is the spectral responsivity of the photosensor. 

In order to formulate the signal current AI ± generated by the spatially varying 
spectral radiant power AP. of the target cone, it is convenient to express the cone 

A 5 i 

radius in terms of the camera instantaneous field of view /3 (see fig. 5(b)) so that 
R = £ /3L, where c is a constant and L is the distance from the target to the camera. 

It is furthermore convenient to assume the two cone areas A + and A_ (as given by 

eqs. (8) and (10), respectively) to be circular in shape with diameters D and D_, 

n 2 + 

respectively, so that A ± = A D ± . Circular symmetry permits the use of the one- 
dimensional Hankel transform instead of the two-dimensional Fourier transform, and 
leads to a significant reduction of computational requirements. This assumption is jus- 
tifiable because most of the radiant power contained in the spatial frequency distribution 
of the cone radiance is generally clustered around the lower spatial frequencies. The 
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assumption is also consistent with the earlier approximation of cone radiance variations 
by two average values . 

The spatial frequency distribution of the target spectral radiant power (in image 
space) may then be given by the Hankel transform of a circular area with radius b ± v a 
as 


ap x( b ±V) = 


= AP 




2n ^ Jq(vs)v dv 


± n2ir pb ± 


= AP 


r ** r u ± v a 

Jo Jo 


dv du> 


where 


„ _ 2?r „ r» _ 27r/ D l' \ r> 

v a A R a n ~ X Ut, ) R a ~ xl 2 r 


2J i( b ± V) 

b v s 

X 


(15) 


and 


s = 


sin O 


k . xffi-U 

i 


± 0L 2R C C lLp 


(16a) 


Upon proper substitution of equations (8b) and (10b), equation (16a) becomes 



(<f>i(e,i, g) > 0 O ) 


b = c 


27 t cos a 


£ COS ej + i(tan ? sh - C sh - tan ^ + £)cos e Q 


“1 1/2 


(16b) 


(<f>i(e,L, g) < 0 Q ) (16c) 


The signal current generated by the photosensor now may be expressed as 

aI ±( u ’ v a; b ±' v ) = J 0 J 0 R x ap x( b ± v a s ) p ( v a s ) L < u > s > J 0 (vs) s ds dx (17) 

where P(v a s) and L(u,s) are the spatial frequency response of the photosensor 
aperture and lens, respectively. 
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The spatial frequency response of the circular photosensor aperture with 


radius v a is 


’M = 


’ c j ° <vs)v 


dv 


_ 2J l( v a s ) 


nV s 

*^0 ^0 


v dv do) 


v a s 


( 18 ) 


The spatial frequency response of the lens is (ref. 10) 

- L(u,s) = cos^| us 2 j|yJ 1 (us) + | sin 2y[j 1 (us) - J 3 (us)| 

- | sin 4y(j 3 (us) - J 5 (us)| + . . ^ 

" A s Sin (^ us2 )f sin y [? 0 (us) " J 2 (us 3 

- | sin 3y[j 2 (us) - J 4 (us)] + | sin 5y£j 4 (us) - J g (us)j - . . A ( 19a ) 


where y = cos 1 . As the photosensor aperture approaches focus (u = 0), the lens 

response becomes 


L(0,s) = — (2y - sin 2y) 


(19b) 


( 20 ) 


The dimensionless variable u is a measure of defocus and is given by (fig. 5(b)) 
u = ^ AH tan fi sin 0 ~ JL AH 
where AH = |ih - # a |. 

Since most of the functions inside the Hankel transform of equation (17) are depen- 
dent on wavelength, the spatial frequency response of the camera cannot, in a rigorous 
formulation, be separated from its spectral response. However, such a rigorous formu- 
lation would lead to unnecessarily difficult computations for the present purpose. It is 
sufficient here to evaluate the Hankel transform only at a single wavelength; namely, the 
surface radiance and camera response weighted average wavelength A 


^ \ S A,i T X,i( t o) p X,i R X,i T X,s,i 
_i 

Yj S A,i T A,i( t o) p A,i R A,i T A,s,i 


( 21 ) 
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Separating spectral and spatial integration in equation (17) and averaging over the 
circular image area of radius b ± v a , this equation may be approximated by (ref. 9) 


AI 


=( u,v a’ b ±) - I 0 

pOO 

^0 


AP x, ± R x dx 


AP X, ± R X dX< 


rb ±Va f H l M L(u,s) J 0 (us) s ds v dv 

Jo Jo Va s v a s 0 


x=x 


I b± 2 v a 2 


f 

2Ji(b ± v a s) 

2 

2Ji(v a s) 

L(u,s) s d^ 

n( v \ 

Jo 

b±v a s 


v a s 

x=x l 2l W 


( 22 ) 

b 2 v a 2 

When b, < 1, then the term — located in front of the spatial frequency integral is 


to be replaced by 


Noise 

In addition to converting the radiation falling on the aperture into an electrical video 
signal and amplifying it, the photosensor and preamplifier also generate random electri- 
cal noise. The root-mean-square magnitude of this noise is designated herein as I g n . 

In the single-line scan mode of the camera, the effective value of this noise can be 
reduced to the value 


*e,n^ n) = 


A e,n 

In 


(23) 


where n is the number of line scans. 

After the electrical signal has been amplified, it is sampled and quantized for 
digital transmission, which causes a so-called "quantization error." For m (=6) bit 
binary encoding, the root-mean-square noise value Iq n (m) of the quantization error is 
(ref. 11) 


l q,n^ = 


*max ~ *min 


\fl2 2 


m 


*max " ^min 
222 


(24) 


where I max and I m i n are the maximum and minimum currents encompassed by the 
camera dynamic range or gain setting. 

Hence, the total root-mean-square value of camera noise may be expressed as 



(25) 
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Image Quality Criteria 


Based on the foregoing results, it is now possible to formulate the three desired 
image quality criteria; namely, scene radiance distribution, minimum detectable albedo 
difference, and minimum detectable cone diameter and slope. Since the spectral depen- 
dence of the illumination scattering function is not important here, it should generally be 
sufficient to evaluate 0^(e, i,g) only at the wavelength X as given by equation (21). 

Estimates of scene radiance distribution .- Important parameters of the scene 
brightness distribution which are to be estimated are the minimum, mean, and maximum 
radiance which can occur as a function of the illumination scattering function, a uniform 
albedo being assumed. Since it is always possible to find surface conditions which gen- 
erate shadows, the minimum signal is defined as (I m i n ) = 0. The mean video signal is 
taken to be that of a flat surface, that is, 

Omean) = *o = ^ ^*X,o^X ^ 

and the maximum video signal is estimated to be that of a surface slope oriented to yield 
the highest possible surface reflectance for a given viewing geometry. The estimated 
maximum level may, therefore, be expressed as 

< I max) = J” p x,ma X R x« < 27 > 

where, from equation (13), 

P X,max = 16^ D d S X T x( l o) p X T X,s^X,max^ ,t,g ) 

Consequently, it is necessary to find the maximum possible value of ^(e,i,g) as a 
function of surface slope a and orientation A brute force method would be to com- 
pute g) for a wide range of a and £, and then to select the maximum value. 

Minimum detectable albedo difference .- Three assumptions are made: (1) the sur- 
face is assumed to be level so that equation (14) can be used to calculate the camera sig- 
nal, (2) the surface albedo is assumed to be constant over the narrow wavelength range of 
the camera spectral filters so that the albedo term can be moved outside the wavelength 
integral, and (3) the albedo difference must be three times the root-mean-square noise 
value in order to be detected. Hence, solving the ratio of signal (eq. (14)) to noise 
(eq. (25)) equal to 3 for p^, the minimum detectable albedo difference may be 

formulated as 
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(28) 


«In 

”^ D /'#'x,oj’” s x T x( 1 o) T x, s R x dx 

Minimum detectable cone diameter and slope .- The signal level of the brighter than 
background area S + and darker than background area S_ may be given from equa- 
tion (22) in the form 



(29) 

The total average signal current excursion is AI = AI + + AI_. 

The accuracy to which this signal can be measured depends not only on the system 
noise but also on the size of the image area. This is simply because the larger the 
image area the more picture elements (or samples of the signal) it contains. According 
to the central limit theorem, the mean of a sample of size n has a standard deviation 
(y/Vn if o is the population deviation. Here a is the total root-mean-square camera 
noise current I^., and the number of independent samples n is the ratio of image area 
over picture- element area, that is, 

„ _ JL 2 ^ h 2) v a 2 .2^2 

n = n b + + b J^2 = b + + »- 

a 

Consequently, the ratio of average signal to root-mean-square noise of a cone 
image may be formulated as 

/ S \ _ AI + + AI - 
\N/" In 

In order to determine favorable lighting and viewing geometries, it should often be 

sufficient to evaluate equation (30) for a fixed cone slope a and normalized cone 
2R C 

diameter c = — — . However, if it is desirable to determine the minimum detectable 

pL 

cone slope, the procedure should be to specify c and then to find the slope a which 
yields a ratio of signal to noise S/N of 3. Or if it is desirable to determine the mini- 
mum detectable cone diameter, the procedure should be to specify a and then to find 
the normalized cone diameter c which yields S/N = 3. Both procedures require 
reiterative evaluations of equation (30) until a solution with desired accuracy is found. 
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CONCLUDING REMARKS 


Image quality criteria were defined and mathematically formulated for the predic- 
tion computer program which is to be developed for the Viking lander imaging experi- 
ment. To define specific image quality criteria, it was stipulated that it is generally 
desirable to observe small spatial details and slopes in monospectral (black and white) 
images and to observe small spectral albedo variations in multispectral (color and infra- 
red (IR)) images. In addition, it was stipulated that it is generally also desirable to 
encompass the entire radiance variations of the scene in all images. Because of the 
limited number of available encoding levels and the potentially wide range of scene radi- 
ance, trade-offs may frequently have to be made between resolving small radiance varia- 
tions and encompassing the complete radiance range. 

To aid in determining these trade-offs for the selection of camera gain and offset 
settings prior to the reception of Mars surface imagery data, three points of the statis- 
tical distribution (histogram) of surface radiance were estimated as follows: A minimum 
radiance was taken to be zero, since shadows will exhibit very little radiance for the 
small atmospheric scattering on Mars in the camera silicon photosensor wavelength 
range. A mean radiance was taken to be that of a level surface and a maximum radi- 
ance was taken to be that of a surface area with a slope and slope orientation which pro- 
duces the highest possible value of the surface illumination scattering function. 

The capability to resolve spatial details and slopes was formulated as the minimum 
detectable diameter and slope, respectively, of a right-circular cone. For steep angles, 
yielding high contrast, this target simulates small detail, and for shallow angles, yielding 
low contrast, this target simulates small slopes. 

The capability to resolve spectral variations was formulated as the minimum 
detectable albedo variation, a level terrain being assumed for computing the lighting and 
viewing geometry dependent surface illumination scattering function. 

An image quality prediction computer program based on the criteria formulated in 
this paper must be carefully compared with images obtained with the Viking lander cam- 
era before being certified for use as operational software. Such comparison may lead to 
changes and additions to these criteria which could improve the value of the operational 
software. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 16, 1973. 
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APPENDIX 


CAMERA CHARACTERISTICS 

Camera characteristics are summarized here to permit preliminary evaluations 
of image quality parameters. Most data are taken from the VLIS Image Quality Analysis 
Report (ref. 12). These data are preliminary and subject to change upon evaluation 
of a developed camera. Some of these data may also have to be modified during the 
mission if camera calibration and engineering data reveal a change in performance 
characteristics. 

Pertinent characteristics of the four camera imaging modes are given in table I. 
The objective lens aperture diameter is = 0.95 cm and focal length is f = 5.25 cm. 
The total average transmissivity of window and optics is = 0.7. 

The electronic noise current is approximately 10 - ^ ampere / /hertz . It should be 
noted, however, that the noise current varies with the amount of radiant power incident on 
the photosensor if the incident radiation is high, and yields a significant increase in noise 
current in the survey mode under some imaging conditions. 

Camera gains and offsets are illustrated in figure 6. The actual gains may be 
slightly different for each one of the photosensors and preamplifiers. These gain differ- 
ences will have to be included in the final prediction program. 

Figure 7 presents a nominal responsivity curve for the silicon photosensor. This 
responsivity may vary significantly, especially in the near infrared (IR), as a result of 
long-term nuclear radiation exposure from the lander radioactive thermonuclear gener- 
ators, and as a function of the photosensor array temperature. These factors must be 
included in the final prediction program. 

Figure 8 presents smoothed transmittance curves for the spectral filters. Final 
curves must be used when they become available. 
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Figure 3.- Cone geometry. 
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Figure 4.- Top view of cone geometry. 
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Figure 8.- Averaged transmittance of spectral filters. 
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